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■ Abstract: We apply the Asymptotic Iteration Method to obtain the bound-state energy spectrum 
' for the d-dimensional Klein-Gordon equation with scalar S(r) and vector potentials V(r). When 
\ S(r) and V(r) are both Coulombic, we obtain all the exact solutions; when the potentials are both 

of Kratzer type, we obtain all the exact solutions for S(r) — V(r); if S(r) > V(r) we obtain exact 

■ solutions under certain constraints on the potential parameters: in this case, a possible general solu- 
te | tion is found in terms of a monic polynomial, whose coefficients form a set of elementary symmetric 
O |. polynomials. 
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I. INTRODUCTION 

The search for exact solutions of wave equations, whether non-relativistic or relativistic, has been an important 
research area since the birth of quantum mechanics. Recently the Asymptotic Iteration Method (AIM) has received 
\ much attention as a method for solving the Schrodinger equation 17 1, both analytically and approximately. It 



o 

has been applied to a large number of physically interesting potentials and has often yielded highly- accurate results. 
Very recently, AIM has been used to study the bound states of the Klein-Gordon and Dirac equations for a number 
of special potentials [I3|-[2^|. In the case of the Klein-Gordon equation, AIM was used to study the bound-states in 
the case of equal vector and scalar potential. To our knowledge, the method has never been used to study the bound- 
states of the Klein-Gordon equation in the case of unequal vector and scalar potentials. In the present paper, we have 
adapted the method to treat problems where the unequal vector and scalar potentials are of Coulomb or Kratzer type. 
. rH , Our main goal is to investigate the exact solutions, whether the scalar and vector potentials, S(r) and V(r), are equal 
or not. In the case of equal vector and scalar potentials, it is known that the Klein-Gordon equation usually reduces 
to a Schrodinger- type equation, which can be studied, for example, by use of the Nikiforon-Uvarov method |30l|. or 
by transforming the equation into a classical hypergeometric differential equation with known solution [3ll]-[32j|. It is 
known [33|-[35| that when S(r) > V(r), bound-state solutions exist. Usually the case when the scalar potential is 
equal to the vector potential is considered separately. The advantage of AIM is that it allows a unified approach that 
can be used to study the bound-state solutions with equal or unequal scalar and vector potentials. 

The Klein Gordon equation has received considerable attention in the literature [35|-[45j]. The present paper is 
organized as follows. Section II is devoted to a brief introduction to the Klein-Gordon equation in (i-dimensions. In 
section III we summarize the Asymptotic Iteration Method (AIM) 1]. In Section IV we use AIM to find detailed 
analytic solutions to the Klein-Gordon equation with Coulomb potentials in (i-dimensions. The lowest even-parity 
solution for the Coulomb problem in one dimension is paradoxical and is still under study (46|-[50T|. Our solutions in 
section IV for d — 1 confirm the exact solutions of the Klein Gordon equation with mixed vector and scalar Coulomb 
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potentials on the half-line obtained earlier by de Castro [48| . In Section V, we use AIM to derive exact solutions with 
equal and unequal scalar and vector Kratzer-type potentials. Our conclusions are presented in section VI. 



II. THE KLEIN-GORDON EQUATION IN d DIMENSIONS 



The d-dimensional Klein-Gordon equation for a particle of mass M with radially symmetric Lorentz vector and Lorcntz 
scalar potentials, V(r) and S(r), r — ||r||, is given (in atomic units H = c = 1) [33j, |3J] by 

{~A d + [M + S(r)] 2 }V(r) = [E - V{r)f^{v), (1) 

where E denotes the energy and is the d-dimensional Laplacian. Transforming to the d dimensional spherical 
coordinates (r, 9i . . . 6d-i), the variables can be separated using 

*(r) = R(r)Yi d _ 1 ,... lh (8i...6 d - 1 ) (2) 

where R(r) is a radial function, and (0i . . . 9<i-i) is a normalized hyper-spherical harmonic with eigenvalue 

1(1 + d — 2), I = 0,1,2,.... Thus, we obtain the radial equation of Klein-Gordon equation in d dimensions by 
substituting Eq.© into Eq.(JTJ) 

- R"(r) - R\r) + | ^ + g~ 2) + [M + 5(r)] 2 - [E — ^(r)] 2 |i?(r) = 0. (3) 

Writing R as R(r) = r _ (- D_1 )/ 2 M(r) gives the radial equation 

- + { (fc '4 r ( 2 fc ~ 3) + [W + S(r)] 2 -[E- V{r)Y] }«(r) = 0. (4) 

where k = d + 21 and u(r) is the reduced radial wave function satisfying u(0) = 0. 



III. BRIEF INTRODUCTION OF THE SOLUTION METHOD 

The asymptotic iteration method was introduced [l[ to obtain exact and approximate solutions of eigenvalue equations 
Q ■ The first step in applying this method to solve Schrodinger-type equations is to transform these equations, with 
the aid of appropriate asymptotic forms, to second-order homogeneous linear differential equations of the general form 

y" = Mr)y' + s (r)y, (5) 

for which Ao(r) and 50(7") are functions in C°°. Here primes denotes the derivatives with respect to r. A key feature 
of AIM is to note the invariant structure of the right-hand side of (J5|) under further differentiation. Indeed, if we 
differentiate (© with respect to r, we find that 



where 



Meanwhile the second derivative of yields 



for which 



y'" = \i(r)y' + Sl (r)y (6) 



Ai — Aq + so + Aq 
s i = s 'o + s oAo- 



y {4) = X 2 (r)y' + s 2 (r)y (7) 



A 2 = X[ + si + A Ai 
S2 = s[ + SqXi. 



Thus, for the (n — l) th and (n) th derivatives of ([5]), n = 1, 2, . . . , we have 

y (n+1) = \ n -i{r)y' + s n -i(r)y 



(8) 

2/("+ 2 ) = \ n {r)y> + s n (r)y 
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respectively, where 



A n — A n _j + S n _i + AoA„_i 
s n = s n -l + s oA n -l- 

In an earlier paper |T[ we proved the principal theorem of the asymptotic iteration method (AIM) , namely 
Theorem 1: Given Xq and sq in C°° , the differential equation f5|) has the general solution 



(9) 



y(r) = exp I - J j-^- 



dt 



C 2 + C! I exp I / (A + ^^)dr I dt 



(10) 



if for some n > 



S n — A n s n _i — A„_is n — 0. 



(11) 



Recently, it has been shown [5l[ that the termination condition (fTTj) is necessary and sufficient for the differential 
equation ([5]) to have polynomial-type solutions. In the next section, we shall apply this method to obtain the exact 
solutions for the relativistic d-dimension Klein Gordon equation with Coulomb potentials. The results obtained in the 
next section are known since the Klein-Gordon equation, in this particular case, reduces into an hypergeometric-type 
differential equation which has known solutions in terms of confluent hypergeometric function. However, the point is 
to unify the technique of using AIM to obtain analytic solutions, by first studying the known Coulomb case, and then 
going on to the more-general Kratzer case. 



IV. THE KLEIN-GORDON EQUATION WITH COULOMB POTENTIALS IN d-DIMENSIONS 

In this section, we use AIM to study the d-dimensional Klein-Gordon equation with vector V(r) and scalar S(r) 
Coulomb potentials [H3]- , namely 

V(r) = ~ S(r) = - S -. (12) 
r r 

In this case, the Klein- Gordon equation (|4|) reads 

- u"(r) + ((*-!)(*- 3) + 4( 8 W) _ 2(Ms + EvU = _ 

In order to solve this Schrodinger-like differential equation using AIM, as mentioned earlier, the first step is to 
transform (fT3"|) into the standard form ([5]). We note that the differential equation (fT3"|) has one regular singular 
point at r = and an irregular singular point at r — oo. The asymptotic solution of (|13|) as r — > oo is given by 
u(r) w e -VM 2 -£ 2 r^ meanwhile the indicial equation of IT3|) at the regular singular point r = yields 

c 2 - c - s 2 + v 2 - -(k - l)(jfe - 3) = 0. 

Thus the exact solution of l[T3|) may assume the form 

u(r) = r c e br f(r) (14) 

where b = —\JM 2 — E 2 and 



c=\ + \J(\-l) 2 + s 2 -v 2 , (15) 

because only the negative root yields a regular wave function at r = 0. Substituting (|14p in (|13[) . we immediately 
obtain 

-2(Ms + Ev) - 2cb 



f(r) . -2 (£ +6 ) fir) + t -W + m-o* /(r) . (16) 
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This is a confluent hypergeometric differential equation which has a regular singular point at r = and an irregular 
singularity at r = oo. The application of AIM is now initiated with Ao = — +6) and sq — ^ 2 ( Ms +^ Ev )^ 2cb ; an( j 
then terminated with the condition (fTTjl . We find 



<5„ = iff [I ( Ms + Ev + cb + kb) = 



k=0 

which in turn yields 

Ms + Ev + cb = -nb, n = 0,1,2,... 
Further, use of (f]~5|) yields the eigenvalue equation 

Ms + Ev 



= n + - + \ <- - l) 2 + s 2 - v 2 . 

VA/ 2 - £ 2 2 V V 2 7 



Upon solving ([19]) for I?, we obtain 



E = M 



sv (3 



(3 2 + v 2 (3 2 + V 2 



\J(3 2 + v 2 - s 2 



where (3 — n + i 



> ., , ., lj + s 2 — v 2 . The un-normalized wave function can be found using (|10p . namely 

s n -i(r) 



/ n (r) = cxp 

which, after some algebraic computation, yields 



A n -i(r) 



-dr 



/„(r) = (-l) n (c + n) n (2c) Bl F a (-n; 2c; 2r v / M 2_ B 2 ) , n = 
up to a multiplicative constant. Here i-Fi denotes a confluent hypergeometric function [5 



,1,2,. 



a a(a + 1) a; 2 

iFi(a;6;a; ^ 1 + r + wiy^ 



(«)n g 
(&)n n! 



E 

fc=0 



(a)fc z fc 
(ft)*, k\ ' 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 



and (a)k — T(a + k)/T(a). Using equations (TT4"]) and (f2"Tj) . we find that the exact solutions of Klein-Gordon equation 
with Coulomb potential in arbitrary dimension d is given by 

u n (r) = a i r5+V(f-i) 2 +- 2 -- 2 e -VM^isV ij p i x + ^(fc_ 2 ) 2 +4(.s 2 -t; 2 ); 2r V / M 2 Te 2 ) (23) 

The normalization constant C„ can be computed by means of J °° \R(r)\ 2 r D ~ l dr = 1, where R(r) = r - ^ -1 " 2 ^?-), 
which, in turn, requires the computation of the definite integral 



Inm(a) = p a e 13 \F\ ( n\ a\ p)\F\ ( to; a; p)dp 



(24) 



when n = to. Since there is some confusion in the notation used in Refs.([53|-[53j|) for the computation of the 
normalization constant C„ in (|22| , we shall present here a detailed computation of the definite integral (|24| . 



Lemma 1: For a > —1, we have 



p Q e p iFi(-n;a; j o)iFi(-TO;a;p)rfp = < 



r(a+n-l) 
(a + 2n)lfM 
[r(o)l 



j/ TO = 71 — 1 

if to = n 



(25) 



±Lmi(n+iy. if m = n + 



l 



I o 



otherwise. 
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Proof: Using the series representation of the Confluent hypergeometric function ([24]) and observing that (— n)& = if 
k > n, we can write ([Ml) 

\ \ - (-rajj (-m) fc /" a+l+k \ - \ - (-njj (-m)fc . 

- g '-^t' + V -l-^a + i+l;a;l) = r,. + l)g H^= , 

as a consequence (5(| of Vandermonde's Theorem 2-F1 (—n, b; c; 1) = (c — b) n /(c) n . The last sum survives for m = 
n — 1, 7i,n + 1. Since we are concerned with the case m — n, we consider 

r(a + 1) <A (-n)i(a + i - l) m 



E 



(a)„ ^ i 



r(a+ 1) f(-rc) n _i(a + l)„_i(-7i)„ (-n) n (a + l)„(-n - l)nl , , n N n!r(a) 
; k — (a + 2n)- 



(a)n { (a) n _i(n — 1)! (a) n n\ J (a)„ 

as required. The other cases follow similarly. O 
Consequently, we have 



/ p a e-P[ 1 F 1 (-n;a;p)] 2 dp=(a + 2n) 1 ^p 



The normalization constant in (|22[) is then provided by: 

C " "Uvm^V ( } (2c) B ■ 

Consequently, the full normalized wave function (|23|) reads 



(26) 



.W - J^S^r ^^(-»; * *^TPj, (27, 

y 2(c + n)n!I (2c) 



. 2 
cases: 



where c = | + \/(| — l) 2 + s 2 — w 2 . Depending on the values of s and v, we may now consider the following three 



• If v = 0, in this case s and M must have the same sign as a result of ([20]) . i.e. s > (since the right-hand side 
of (Til?)) is greater than 0). Thus we have 

( .* V' 2 

E = ±M 1 . (28) 

V (n+i + ^(f-l) 2 + S 2 )V 

• If s = 0, then, again using (fT5|) . i? and u must have the same sign for (-| — l) 2 > v 2 . For attractive Coulomb 



2 

-U / — I \W ' I W 1 1-1 * ,1 1 1 1 '( 1 "7 _ II I 1 <•/ 1 II I / ' tlhlltl \ 1\ 

1/2 



potentials < v < § + 1 — 1, we have the so-called Tr-mesonic atom [5 



E = M\1+ ] . (29) 

( n+ l +v / ( |_l ) 2„„2 ) 2 



where n and I are the radial and angular quantum numbers. 
• If s = v > 0, since \E\ < M, we have 

2v 



2 



E = M \l - -— . (30) 

V (»+!-§) + v v 
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V. THE KLEIN-GORDON EQUATION WITH KRATZER POTENTIALS IN d-DIMENSIONS 

In this section, we consider the scalar and vector potentials with the forms 

S ( r ) = -- + ^> V ( r ) = -- + ^- ( 31 ) 

The Klein-Gordon equation ^ then reads 

„ (-2(M Sl + Ev x ) 2Ms 2 + s 2 1 +2Ev 2 -v 2 + \(k-l){k-3) 2( Vl v 2 - Sl s 2 ) . s 2 2 -v 2 2 \ fJp2 ^ 
— u + < 1 t. 1 5 1 -. — >u={hj — M )u. 

(32) 

This differential equation has two irregular singular points at r — and r — oo. As r — > oo, the differential equation 
u" ~ (M 2 — E 2 )u has a solution given by u » exp(— yP — E 2 r). As r — ► 0, it is clear that the differential equation, 
after replacing z = 1/r and letting z — * oo, has a solution it rs exp(— ys| — vf r _1 ). Thus we may assume the exact 
solution of (l32l) takes the form 



/l2 

V s 2~"2 



Substituting this expression in (|32p . we obtain 

-2(Msi + 2(Ms 2 + £fa 2 ) + - v\ + \{k - l)(k - 3) + 2a6 



g"{r) = -2(-- s +b)g'(r)+l — ^ - '- + 

2{v\v 2 - sis 2 ) - 2a 



9{r), (34) 



where we denote a = —\js\ — v\ and b = —\JM 2 — E 2 . 



A. Equal scalar and vector potentials 

In order to solve (|3"4"|) , we consider first the case of equal scalar and vector potentials, i.e. V\ = s\ = B and v 2 = s 2 = A. 
Then, equation ([34]) reduces to 

»( \ oh 'Mm f -2B(M + E) 2A(M + E) + \(k-l)(k-3) \ 
g (r) = -2bg (r) + + 1 \g[r). (35) 



This equation has a regular singular point at r — with indicial equation given by c 2 — c— 2A(M +E) — \{k— l)(k— 3) 
0, which implies 



c=^ + y(^-l) 2 + 2A(M + i?), (36) 



where we assume c > (because only the negative root yields a bounded solution at r = 0.) With g(r) = r c f(r), 
equation (|35p now reads 



/"(r) = -2 (£ + &)/» + (: 



■2B(M + £7) - 2c6 



/(r)- (37) 



By a direct application of AIM with Ao = f + b and sq = 2B ( m +e) icb ^ ^ e ^ ermma y on condition (jlip implies 
S n = 0, n = 0, 1, 2, 3, . . . if and only if 

J] (B{M + E)+cb + kb)=0,^c + n= ^==== - (38) 

Combining (|36p and (|38p . we obtain the eigenvalue equation 

7m=w = 2n + 1 + ^ (fc ~ 2)2 + 8A(M + ^ (39) 
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which generalizes the results of [HJ for the 3-dimension case. Although we can solve Eq. ((39j) exactly for E, the 
expression for E is rather complicated. We therefore write Eg. ([55]) as 



2By/2M - (3 2 - [3^2n + 1 + \J (k - 2) 2 + 8A(2A/ - /? 2 )^ =0 



(40) 



and expand it as a power series in /3 = \JM — E. Thus we have 

B 



2Bv / 2M + (-2n - 1 - = 2 ) 2 + 16AM)/3 - -^=P 2 + • • • = 0. (41) 

As an approximation, if we take only the term in j3, and solve the energy equation, we obtain the non-relativistic 
energy [57J] 

E' = E-M = -P 2 = 8MB2 — (42) 

(2n + 1 + v/(fc - 2) 2 + 16AM) 2 



The exact solutions f n (r) of (l37|) can be obtained using (fT0|) to yield 



/„(r) - (-l)"(c + n) n (2c) nl F 1 (-n; 2c; 2^M 2 - E 2 r) (43) 



up to a constant. Consequently, the exact solutions (|33|) of ([32)1 . in the case equal scalar and vector potentials, are 
given by 

Un (r) = N n r c e~ VWriWr 1 F 1 (-n;2c; 2ry/M 2 ~E 2 ), n = 0,1,2,... (44) 
where the normalization constant (implied by Lemma 1) now becomes 



(2VM 2 -E 2 ) 2 ^(2c) n 
y 2(c + n)ml (2c) 



where c is given by (|36p . We note that in the case A = 0, equation (f3"5)) yields 



m + £ - ir + £ 2 



2 2 

which is in complete agreement with our results for the Coulomb potential. 

B. Unequal scalar and vector potentials 

In the case of different scalar and vector potentials, we have a ^ and therefore we have to approach equation 
directly. The equation, however, has two irregular singular points, at r = and r = oo, and the standard techniques 
of solving differential equations cannot be applied. AIM, however, has the advantage that it can be applied directly to 
obtain exact solutions under certain conditions on the potential parameters of (|34"|) . Indeed, with Ao(r) = — 2( — -pz + b) 

an d SQ ( r ) = + 2iMs a+ Ev a ) + s' 1 - V l + i(k-l)(k- 3 ) + 2a b + ^-s^-Za ^ termination rdation JJJJ 

for 5 n = 0, where n = 1, 2, . . . is the iteration number, the following conditions 

2nb = -2{Ms 1 + Ev 1 ), 

n(n - 1) = 2{Ms 2 + Ev 2 ) + sj - vj + |(fc - l)(fc - 3) + 2ab (47) 
—2na = 2(v±V2 — S1S2) — 2a. 
Furthermore, the exact solutions under the constraints (|47|) . using the AIM expression (fT0| . are given by 



(r) = r" u(r) = C n r n e~^ M ' 2 - E " r -^^ , n=l,2,.... (48) 
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The normalization constants C n of (|48|) can be computed by the standard identity [58[ (in which K v is the modified 
Bessel function of second kind of order v) 



which implies that 



. c + 1 



•~ c exp[-£r- ^}dr = 2( — ) Jf_c_i(2V-4B), .4>0,£>0. (49) 



C » 2 = 2 ( ) ^ " L(^-^>W-r^.. (50! 



That is to say, the exact solutions for the differential equation (|34[) under the constraints (|4"T|) are given by 



1+2 



un(r) = [2 ( 4 ^-x(2(^ - «f)*(M 2 - E 2 )i)j r « e -^/^-^^P. (5 1) 
Further, we note that the first and the third equations of (|4"T|) yield 



_ t>l-»2-SlS 2 _|_ ^ 



(52) 



This, in turn, inplies 



<s\ - v 2 y/M 2 - E 2 = (Msi + Evi)\jsl -v 2 - (viv 2 - s 1 s 2 )V ' M 2 - E 2 . (53) 
Thus, by mean of the second equation of (|47p . the eigenvalue equation is 
Mai + Evi 1 



VM 2 - E 2 2 V v 2 



Some important remarks are in order: 



(- - l) 2 + 2(M«i + E Vl )Js 2 -v\- 2{ Vl v 2 - Sl s 2 )^M 2 -E 2 + 2{Ms 2 + £u 2 ) 



(54) 



• The exact solutions it„(r), n = 1,2,... are nodeless for all n. Here, n represents the AIM iteration number. 
Thus, for all the values of the potential parameters that satisfy (p)Tj) - (f52")) . the corresponding solution represents 
a ground-state wave function. Note that each operator now acts on a different parameter space. 



• We notice for n — 1, the second equation of (|5ip implies v^v 2 — sis 2 and since s 2 > v 2 , we must have vi > Si 
which in turn implies — si > —v±. Thus, the requirement for bound states S(r) > V(r) is satisfied. It is also 
clear that, this is the case for all n > 1. 



• If s 2 =v 2 = 0, we obtain u Q {r) = C rHV(|- 1 ) 2 + s ?-«'?e-^- £2r and 



Msx + Evt 1 Lk 



V M 2 - ff 2 2 ' V ( 2 1)2 + S? ^' 
which reduces to the ground state solution described for Coulomb potential as mentioned earlier. 
• In the case of a pure scalar potential, i.e. V(r) = 0, we have v\ — v 2 = 0, S(r) = — — + |J > 0, 

u (r) = Cori + ^^^ 1 ^ 2+2MsiS2+2si32VMJ ~ E2+2M32+s ^e~ VM ^ E ^ r ^^' 

and 

M Sl 1 



VM 2 - E 2 2 V 2 
If we further assume that s 2 — 0, we obtain 



(- - I) 2 + 2M Sl s 2 + Sl s 2 \/M 2 -E 2 + 2Ms 2 



1/2 



which again agrees with the result for the Coulomb potential for this particular case. 
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In order to move beyond the ground state solutions, equation (|5Tj) suggests that the exact solutions u(r) of ([34]) take 
the form 



U(r) = r c e -^M^Wr-^l f{r y 



(55) 



Substituting this expression in (|33p . we obtain 



f"(r) 



,c a 



b)f'(r) 



-2(Ms! + Evi) - 2cb 2(Ms 2 + Ev 2 ) + si -v\ + \{k - l)(k - 3) - c 2 + c + 2ab 



r r 

2{v\V 2 — S1S2) — 2a + 2ca 



/(r) 



where again a = — -^/s 2 — v 2 and 6 = — \J M 2 — E 2 . We may now choose c such that 



w i^2 — S1S2 — a + ca = u ^ c = — : + 1. 



which in turn reduces Eq. ((56|) to 

f"(r) = _ 2 (£-^ + 6)/'(r) + (: 



■2(M Sl +Evx)-2cb G\ f 
+ — ) /(r), 



where we denote 



G = 2(Ms 2 + Ev 2 ) + sj -v\ + -{k - l){k - 3) - c 2 + c+2ab. 



(56) 



(57) 



(58) 



(59) 



Although (15811 still has an irregular singular point at r — 0, the direct application of AIM with Ao = — — ^ + 6) 
and sq = ~ 2 ( Msi +- E,,1 )~ 2cb _j_ implies, by means of the termination condition (fTT|) , that 



Msi + Sui + cb = -nb 



c + n 



Ms 1 + Ev x 



VM 2 - E 2 

where for n — 0, 1, 2, . . . the following constraints on G must hold 



n = 0,1,2,..., 



(60) 



So 



= G = 

= G 2 - 2cG - 4a6 = 

= G 3 - 2(3c + 1)G 2 + 4(2c 2 + c - 46a)G + 16a6(2c + 1) = 

= G 4 - 4(3c + 2)G 3 + 4(1 lc 2 + 3 + 13c - 10ao)G 2 + (192afe - 72c 2 + 240cao - 24c - 48c 3 )G 
- 1446a + 144a 2 6 2 - 432ca6 - 288c 2 a6 = 



and so on, for higher iteration numbers. 



C. The case Go = 

In the case G = 0, we have 



c - X - + \J (| - l) 2 + 2^s 2 ~v 2 ^jM 2 -E 2 + 2(Ms 2 + Ev 2 ) + s 2 - v 2 (61) 

as a root of (l59|) and, again, c > because only the negative root yields a regular wave function at r = 0. On the 
other hand, equation ((60|) implies 

Msi + Evi 

c = ■ , (62) 

VM 2 - E 2 
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where now ([57)) gives 



'4 - v 2 yjM 2 -E 2 = (Msi + Evx^sl -v\- (viv 2 - s\s 2 )^/ M 2 — E 2 . (63) 
Therefore, the ground-state energy equation in d-dimensions is given by 

m j^ Ev e 2 = l + \Af~ 1)2 + 2(Msi + - 5 - 2{viv2 - s ^ m2 - e2 + 2 ( Ms2 + Ev ^ + s ? - «?> 

(64) 

which in complete agreement with equation ([54]) . 

D. The case Q x = 



In this case, the constraint C*i := G 2 — 2cG — Aab = implies G = c ± Vc 2 + 4a6, which implies by ([55)) for G > 0, 
that 

y/c 2 + 4a6 = 2(Ms 2 + Ev 2 ) + s\ - vj + -{k - l)(fc - 3) - c 2 + 2a6. (65) 
In order to simplify the notation, we denote 

fx = 2(Ms 2 + Ev 2 ) + s\ - v\ + j{k - l)(k - 3) + 2ab (66) 



and solve (|6"5"|) for c > to obtain 



c= iy / 2 + 4^ + 2 v / 47T+l + 16a&. (67) 



Meanwhile, ([57]) and J60]) imply, for ab = y 'M 2 - E 2 ^Js\ - v 2 , that 

a& = - (Msi + Ev x )\l s 2 2 - v\ 
The eigenenergy is then given, from (|67|. by the equation 



ab = ~(M Sl + Ev 1 )^s 2 -v 2 -'^{v 1 v 2 - Sl s 2 )^M 2 -E 2 . (68) 



Msi+Evi 1 



y/M 2 - E 2 2 



2 + 4/i + + 1 + 16a6 - 1, (69) 



where fi and a& are given by ([66]) and ([68]) . respectively. Further, the exact solution of ([56]) and G = c + y/c 2 + 4ab 
becomes 

/iW=r-| (70) 

up to a constant. Thus 



ui(r) = d(r - | )r c e -yiCTr-^P ; (n) 
where the normalization constant C\ can be computed by means of (|49[) . 

E. The case 2 = 

In this case, we have for Ms x + Ev x + cb = -2b that Q 2 := G 3 - 2(3c + 1)G 2 + 4(2c 2 + c - 46a)G + 16a6(2c + 1) = 0, 
where G is given by (|59[) . The eigenvalue equation is then given by the root of this equation along with G given by 
([59]) and 

ab = i(M Sl + Evx^sl-vl- i(«!«2 - s lS2 )VM 2 -E 2 , (72) 

in a similar fashion to the previous case. After some algebraic computations, we obtain for the exact solution that 
G-4c-2 1 r—— ■ — — — w G-4c-2 1 



/a(r) = (r + — y/lQab + 8c + 4 + AcG - G 2 )(r — — yjl&ab + 8c + 4 + 4cG - G 2 ). (73) 

4o 4o 46 4o 
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F. The general case Q n := 



The structure of the wave functions for the cases Go = 0, Q\ = 0, and Q 2 — 0, etc., allows us to formulate a possible 
general solution of the Klein-Gordon equation in the case of unequal scalar and vector Kratzer potentials for arbitrary 
n and d. For this purpose, we assume that the general form of the exact solution of the differential equation (|59|) 
takes the form of a monic polynomial 



Mr) = H(r-a k ) = 



k=l 



k=a 



(74) 



whose coefficients {a k }^ =0 are the elementary symmetric polynomials [5£ 

T-n = 1, 



O-n-1 



«n-2 



«n-3 



(75) 



a k 



a = 



On substituting (fT4"| into the differential equation 
set of linear equations: 



l<k<n 

J2 cr i cr J) 
l<i<j'<n 

l<i<j<k<n 

E ^1 • ■ ■ a j k 

l<Ji<72<— <jfc<" 
n 

(-1)" n ^ 
fc=i 

we find that the coefficients {fifc}J? =0 must satisfy the following 



( Ga + 2aa x = 0, 

4aa 2 + (G - 2c)a x + 2nba = 0, 

[n(n — 1) + 2cn — G] a n = 2ba n -\, 

k 2o(A; + 3)a k+3 + (G - 2c(k + 2) - (k + l)(k + 2))a k+2 + (F - 2b{k + l))a k+1 = , for < k < n - 3. 
In order to understand the application of these formulae, we consider the case of n = 3: equation (76) implies 



(76) 



f a 2 = -i- a (G-2c 



Ylab- 



Ol, 



6+6c-G a2 ' 

6aas + (G — 4c — 2)a 2 + 46ai = , since k 
The solution of this system, for a x ^ 0, gives the condition Q 3 on G that reads 



(77) 



0. 



G 4 - 4(3c + 2)G 3 + 4(llc 2 + 3 + 13c - Wab)G 2 + (192ab - 72c 2 + 240ca6 - 24c - 48c 3 )G 
- 1446a + 144a 2 6 2 - 432cafe - 288c 2 a6 = 



That is to say, the same condition that we obtain if we apply AIM to (|58[) with three iterations, i.e n = 3. Equations 
(|76|) are the full set of restrictions on the parameters of the Kratzer potentials (|3T|) in the case of unequal scalar and 
vector parts. 



VI. CONCLUSION 



In this paper, the Klein Gordon equation in arbitrary dimension has been solved exactly for the bound states corre- 
sponding to Coulomb or Kratzer scalar S(r) and vector V(r) potentials. When both potentials are Coulombic, equal 
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or not, we find all the analytic solutions. When both potentials are of Kratzer type, we find the exact solutions when 
S(r) = V(r); when the potentials are unequal and the scalar potential dominates, we generate exact solutions under 
certain specific conditions on the potential parameters. Furthermore, a general solution is found in terms of a monic 
polynomial whose coefficients form a set of elementary symmetric polynomials. Our method of solution is based on 
the recently-introduced Asymptotic Iteration Method. This approach has the advantage of simplicity in the exact 
cases, and flexibility, leading to approximations, when exact solutions are not attainable. 
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